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Abstract 

We consider the Krall-Sheffer class of admissible, partial differential operators in the plane. 
We concentrate on algebraic structures, such as the role of commuting operators and symmetries. 
For the polynomial eigenfunctions, we give explicit forms of the the 3— level recurrence relations 
and differential raising operators, which are shown to satisfy unusual commutation relations. 
We present new generating functions for two of the cases. 
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1 Introduction 

Orthogonal polynomials appear in many places in mathematics, physics and engineering. In math- 
ematical physics, orthogonal polynomials arise in the context of eigenfunctions of linear, second 
order differential operators. Bochner (see [T], pl50) classified those operators which have exactly 
one polynomial eigenfunction of each degree. The classical, one-dimensional examples possess im- 
portant structures which enable their explicit construction. Among these are the 3— point recurrence 
relations and the differential ladder operators, which allow the polynomials to be built in a recursive 
manner. Starting from Rodrigues formulae and generating functions, these structures can themselves 
be built in a systematic way. 

The classification of orthogonal polynomials in the plane and higher dimensional spaces is much 
more complicated than the case of the real line. However, requiring the polynomials to be eigen- 
functions of linear differential operators makes the classification problem more feasible. An early 
and systematic treatment of such polynomials in two variables, satisfying a second order eigenvalue 
problem, was given by Krall and Sheffer |S]. They showed that, up to affine transformation, there are 
9 distinct cases. Three of these just correspond to products of classical 1— dimensional polynomials, 
but the remaining six are non-trivial extensions. 

To generalise Bochner's result for the classical orthogonal polynomials, Krall and Sheffer [8] 
introduced the concept of an admissible operator. In two dimensions we have to expect the existence 
of N + 1 polynomial eigenfunctions of degree N. 
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Definition 1.1 (Admissible Partial Differential Operator) The differential operator L is said 
to be admissible if, for the eigenvalue problem 



Lip = Xip, 

there exists an infinite sequence of eigenvalues {Aat}^ =0; such that for each A = An there are no 
non-zero polynomial solutions of degree d < N and there are exactly N + 1 linearly independent 
polynomial solutions of degree N . 

It is easily seen that a second order admissible differential operator must have the form 



Lip := (ax 2 + d\X + e x y + f\)ip xx + (2axy + d 2 x + e 2 y + h)<p xy 



+(ay 2 + d 3 x + e 3 y + f 3 )<p yy + {]3x + Ki)ip x + (Py + K 2 )ip y . 

The number of parameters is reduced further by making affinc transformations on the x — y space, 
giving 9 canonical forms (see Section [5]) . The degree TV polynomial cigcnfunctions form an N + 1 
dimensional vector space Vm and we may choose a basis of monic polynomials: 

Pm.n = x m y n + lower order terms, for m + n = N. (2) 

The eigenvalue for P m ^ n is easily calculated from the leading order term to be 

Kn+n = (m + n)((m + n - I) a + 0). (3) 

We see that when j3 is not a negative integer, A mi + ni 7^ A„ l2+ „ 2 for mi + ri\ 7^ ra 2 + n 2 . 

This family of polynomial eigenfunctions can be arranged into a triangular array (see Figure 
[T]), with each horizontal level representing a particular Vjv- In place of the three point recurrence 
relations for the one dimensional case, Krall and Sheffer showed that in two-dimensions there are 
two three level recurrence relations of the form 

m+n m+n— 1 

P'm+l.n — xP 7n n -\- ^ O^i^m^nPi.m+n—i ~\~ ^ ^ bi,m,nPi,m+n—l—ii 
i=Q i=0 

(4) 

m+n m+n — 1 

Pm,n-\-l yPrn.n ~\~ ^ ^ ^i,m,nRm-\-n—i,i ^ ^ di rn n P 7n j rn _ i_ ^ ^ , 

i=0 i=0 

but gave no specific forms, apart from the case of direct products of classical one dimensional 
polynomials. Such 3 level recurrence relations have recently been discussed more generally in [3J|2]. 
A detailed discussion of orthogonal polynomials in two variables can be found in 

The coefficients of the second order terms of the operator ([T]) can be considered as the components 
of an upper index metric tensor, with L being a deformation of the corresponding Laplace-Beltrami 
operator L^. In |12[ [6] it was observed that each of the 9 metrics was either flat or constant curvature 
and that the additional first order terms correspond to "trivial" magnetic field, which can therefore 
be gauged to a purely potential term. They also showed that each of these cases (considered as a 
quantum system) was super-integrable and presented two commuting operators for each case. 

In this paper we first discuss the symmetries of these Laplace-Beltrami operators and use these 
to group together the 9 canonical forms into 4 classes (see Section [2]) . Within each class the metrics 
are related through (non-affine) coordinate transformations. Each metric is invariant with respect 
to an involutive transformation, corresponding to an involutive automorphism of the algebra of 
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symmetries. These involutions are very useful in Section [31 where we discuss super-integrability. 
We briefly review separable potentials, defining L = Lb + U(x,y), each of which is characterised by 
the existence of a commuting, second order operator and each of which depends upon two arbitrary 
functions of one variable (as in the classical mechanical case; see [5]). Requiring two such commuting 
operators (for super-integrability) reduces these potentials to simple rational functions depending 
upon (at most) 3 arbitrary parameters. Super-integrability (with linear or quadratic commuting 
operators) is associated with separability in more than one coordinate system (see, for example, [7]). 
It is then possible to gauge transform these super-integrable cases to Krall-Sheffer form (just the 
reverse calculation of [5]), but this sometimes restricts the parameters to a sub-case. 

Our main results are in Sections 2] and [5] In the first of these, we discuss the algebraic structure 
of the polynomial eigenfunctions of the Krall-Sheffer operators. We first emphasise the connection 
between super-integrability and the construction of eigenfunctions. The commuting operators can be 
used to build the two-dimensional polynomial eigenfunctions out of one-dimensional polynomials, 
satisfying a reduced differential operator. In Section HJ we also give explicit formulae for the 3- 
level recurrence relations for all the (non-trivial) Krall-Shefffcr operators. We also present explicit 
formulae for a pair of differential ladder operators, which satisfy R+, x P r ,i,n = Pm+i.n, R+,yPm,n = 
Pm.n+i {raising operators). None of these formulae have previously appeared in the literature. 

In Section [5] we present generating functions for the polynomials of Cases V and VIII (both new) 
and of Case IX. The relationship (discussed in Section [5]) between Cases I and IX, gives us an array 
of eigenfunctions (not all polynomial) of a restricted form of the Case I operator. 

2 The Geometry of Krall-Sheffer Operators 

For an n— dimensional (pscudo-)Ricmannian space, with local coordinates x , • • ■ ,x n and metric <?y, 
the Laplace-Beltrami operator is defined by L^f = g 13 'VjVj/, which has explicit form 




where g is the determinant of the matrix gij . The coefficients of leading order terms in the Laplace- 
Beltrami operator are the coefficients of the inverse metric g 13 . For a metric with isometrics, the 
infinitesimal generators (Killing vectors) are just first order differential operators which commute 
with the Laplace-Beltrami operator (JSJ) - When the space is either flat or constant curvature, it 
possesses the maximal group of isometries, which is of dimension |n(n -I- 1). In this case, Lb is 
proportional to the second order Casimir function of the symmetry algebra (see [5]). The Krall- 
Sheffer metrics fall into this category, with n = 2. 

When n = 2 we have 3 independent Killing vectors for flat and constant curvature metrics. 
For the particular metrics defined by each of the Krall-Sheffer operators, a basis of Killing vectors 
is easy to calculate explicitly. Four Lie algebras occur and the 9 canonical forms can be grouped 
under these, with the members of each group being related through nonlinear transformations of the 
coordinates (see Table [TJ. The algebras so(3) and si (2) generate metrics with constant curvature, 
while the Euclidean and pseudo-Euclidean algebras generate flat metrics. For the purposes of this 
grouping we label the coordinates (xj, j/j), i = 1, . . . , 9. 

2.1 The Rotation Algebra so(3): Cases I and XI 

The abstract algebra, with A, B and C satisfying 

[A,B] = C, [B,C]=A, [C,A]=B, has Casimir function C = A 2 + B 2 + C 2 . 
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Cases 


Symmetry algebra 


I, IX 


so(3) 


II, III 


sl(2) 


IV, VI, VII 


e(2) 


V, VIII 


e(l,l) 



Table 1: The symmetry algebras for the Krall-Sheffer metrics 
The concrete realisation 



A = - xi - yi)d Xl , B = 2 y / y 1 (l -x\- yi) d m , C = (d Xl - d yi ), 

gives the Laplace-Beltrami operator 



L 6 = —C = (xi - x l )di 1 + 2x m d Xl d yi + (yf - yi )d 2 yi + - ((3xi - 1)^ + (3j/i - 



The concrete realisation 



A = y 1 - £9 - 2/g d X9 , B = y 1 - xl - yl d yg , C = j/ 9 a xg - x 9 d Vg , 

gives the Laplace-Beltrami operator 

L b = -C= (x 2 9 - l)8 2 Xo + 2x 9 y 9 d X9 d yo + (y 2 9 - l)d 2 y9 + 2(x 9 d Xo + y 9 d yg ). 

In each case, the involution 

(x,y) = (y,x) gives the automorphism A <H> B, C — >• C, 

which leaves the Lf, invariant. The two realisations are related through the change of coordinates 
xi = xl, yi = yl. 

2.2 The Special Linear Algebra sl(2): Cases II and III 

The abstract algebra, with E, F and H satisfying 

[H, E] = 2E, [H, F] = -2F, [E, F] = H, has Casimir function C = H 2 + 2EF + 2FE. 
The concrete realisation 

r V2 



U = Ax 2 d X2 , E = 2^x^d y2 , F = 4^x^d X2 +2(y 2 -l)J^d y2 , 
gives the Laplace-Beltrami operator 

L b = -^C = x\d 2 X2 + 2x 2 y 2 d X2 d V2 + (y 2 - y 2 )d 2 2 + -(3x 2 d X2 + (3y 2 - l)0 y2 ), 
which is invariant under the involution 

(x 2l y 2 ) = ( — —,2/2 ) , giving the automorphism E <-> F, H — > H. 

V x 2 J 



The concrete realisation 



H = ix 3 d X3 + 2y 3 d ya , E = x 3 d y3 , F = iy 3 d X3 + 4x3 + 3 ^ 3 d ya ; 

X3 



gives the Laplace-Beltrami operator 

1 3 

L b = = x\dl 3 + 2x 3 y 3 d X3 d V3 + (y\ + x 3 )d; 3 + -{x 3 d X3 + y 3 d V3 ), 

which is invariant under the involution 

[x 3 ,y 3 ) = ( / (j/3+ 3 X3)2 , fal + 4^ 3 )z/ 3 \ ^ g . ving the automorphism e -s-> F, H -> H. 

V x 3 x 3 / 

The two realisations are related through the change of coordinates X2 — —64a;3, 7/2 = — 2/3/(4x3). 

2.3 The Euclidean Algebra e(2): Cases IV, VI and VII 

The abstract algebra, with A, B and C satisfying 

[A, B] = 0, [C, A] = B, [B, C] = A, has Casimir function C = A 2 + B 2 . 
The concrete realisation 

A. = ^/x~ld Xi , B = y/y^d yi , C = 2y/x 4 y4 (d Vi - d x J, 
gives the Laplace-Beltrami operator 

L b = C = x 4 d 2 X4 + y 4 dl + ^(d X4 +d y4 ), 
which is invariant under the involution 

(X4, 274) = (2/4j X4), giving the automorphism A -f-> B, C — >• C. 
The concrete realisation 

A = y/xEd X6 , B = d Ve , C= y/x^ (2d V6 - y 6 d Xe ) 
gives the Laplace-Beltrami operator 

L b =C=X 6 d 2 x6 +d 2 y 6 + ±d X6 , 

which is invariant under the involution 

(x6,2/6) = > gi vm g t ne automorphism A -s-> B, C — > C. 

The concrete realisation 

A = <9 X7 , B = a y7 , C = x r d y7 - y 7 d X7 
gives the Laplace-Beltrami operator 

Lb = C = d 2 X7 +d 2 U7 , 

which is invariant under the involution 

{x7,yi) = (2/7, X7), giving the automorphism A <-> B, C — > C. 

These realisations are related through the changes of coordinates: x 4 = xq, y 4 = z/|, and xg = 
xf , 2/6 = 2/7- 
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2.4 The Pseudo-Euclidean Algebra e(l,l): Cases V and VIII 

The abstract algebra, with E, F and H satisfying 

[H, E] = E, [H, F] = -F, [E, F] = 0, has Casimir function C = EF. 
The concrete realisation 

H = 2x 5 d X5 , E = y/x^d V!il F = 2^/x^d X5 + —— d ys , 

V x 5 



gives the Laplace-Beltrami operator 



L b = C = 2x 5 d X5 d V5 + y 5 d 2 + d V5 . 



which is invariant under the involution 



(^5,2/5) = ( —,U5 ) giving the automorphism E -f-> F, H — > H. 



£5 



The concrete realisation 



4 

gives the Laplace-Beltrami operator 



H= I x 6 - ^y'i) d xs -y 8 d ys , E = y s d Xa +2d ys , F = d Xg , 



L b = C = y & d 2 + 2d Xa d Vs , 



which is invariant under the involution 

(x 8 ,y a ) = ^I^L±if^L + Hi j 2x s - giving the automorphism E <H> F, H -> -H. 

These realisations are related through the changes of coordinates: 3:5 = ^ (4x§ — y|) 2 , 2/5 



3 Adding Potentials: Separability and Super-Integrability 

Second order operators, commuting with Lf,, are just symmetric quadratic forms of Killing vectors. 
Suppose K is such an operator. Then we may seek functions U and V, such that 

[L b + U, K + V] = 0, (6) 

which constitutes a coupled system of partial differential equations for U and V. The solution 
depends upon a pair of arbitrary functions, each of one variable, typical of separable systems. 
Requiring that there exist two of these second order commuting operators, Ij = Kj + Vj, strongly 
constrains these (formerly) arbitrary functions, which reduce to rational functions depending upon 
only a finite number of parameters. The coefficients of the second order derivatives in such an 
operator K define a contravariant, rank-two Killing tensor. For brevity, we will refer to these 
operators as Killing tensors in what follows. 

Such calculations are standard, so we only present (brief) details of Case I and then just some 
super-integrable cases for the others. In the calculations it is only necessary to consider one of each 
of the equivalence classes listed in Table [TJ Within each equivalence class we can (approximately) 
halve the number of choices of quadratic Killing tensor, by using the corresponding involutive auto- 
morphism. 
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3.1 Krall-Sheffer Metrics I and IX 

First consider Ii = A 2 + V in ([5]). The bracket of two second order operators would generally 
be third order, but in this case the third and second order terms automatically vanish, since A 2 
commutes with Lt>. For the remaining 3 coefficients to vanish, we require 

(x - l)V x + yV v = 4(1 - x - y)U x , xV x + (y - 1)V V = 0, L b V = A 2 U. 

The intcgrability condition for the first pair is just 

d x (xU x + (y-l)U y + U) = 0, 

leading to the general solution of both U and V depending upon 2 arbitrary functions: 

U(x,y) = U 1 (y) + U2 ^\\ V(x,y)=4U 2 (^). (7) 
y-1 V x J 

The third condition is then automatically satisfied. 

Remark 3.1 (Separation of Variables) Written in terms of the coordinates £ = (y — l)/x, ?/ = 
y, the operator L (with the potential ^7ty) is separable. 

The case of I2 = B 2 + V is obtained from 1\ by using the involution. In this particular case, this 
is a very elementary operation of switching the roles of x and y. 

The case of -Z3 = C 2 + V is calculated independently and leads to another second order PDE for 
U(x,y): (xd x + yd y + 2)(U X - U y ) = 0, giving 



U(x,y) = U 3 [x + y) + -^-, V(x,y) = -AuJ-) , (8) 
x + y \xs 

which is symmetric (in form) under the involution. 
3.1.1 Super-Integrability 

Requiring two of above operators (say I\ and I2) to commute with L reduces the potential to 
a rational function. In this case we have three such commuting operators, L = Lb + U, Ii = 
A 2 + Vi, I2 = B 2 + V 2 , h = C 2 + V 3 , where 

j j = fc l , ^2 ^3 v j kl x , 4fc 2 (y-l) 

(1 - x - y) + x + y ' Vl (*+y-V - ' 

_ Ak lV 4k 3 (x-l) _ 4k2y 4k3X 
V2 — 7 ; 7T H , v 3 — 



(x + y-1) y 



These operators are not independent, satisfying AL + 1\ + 12 + 13 = 4fci , which reduces to 4L^ + A 2 + 
B 2 + C 2 = 0, when the parameters are zero. If we extend the involution to act on the parameters 
as fci hi, &2 k 3 , then 

L^L, h*+h, h->h- 
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3.1.2 The Transformation to Case IX 

We have a transformation between the Laplace-Beltrami operators of cases I and IX. The potential 
([7]) is transformed to 



X 2 J 



U 2 [^) f.* 



U(x,y) = U 1 (y 2 ) + -L, V(x,y) = 4U 2 

y - 1 



ir-i 



The potential ((HJ is transformed to 

U(x,y) = U 3 (x 2 +y 2 ) + \ 1 , ^(x,y) = -4(7 4 % 

x + y \x z 

In particular, the potentials © give us the corresponding supcr-intcgrable case for Case IX. 
3.1.3 The Krall-Sheffer Gauge 

Starting with the operator L = + U , with U defined by ([3]), the gauge transformation 

L i-^ G~ X LG, with G = (1 - x - y) a x b y c 

leads to k\ = \a(2a — 1), k 2 = \b(2b — 1), fc 3 = ic(2c — 1) if we require the potential to reduce to a 
constant. The operator then reduces to the Krall-Sheffer form with 

o=-(208 + «i + «2)-l), 6 = -7(2«i + l), c=-i(2rt 2 + l)- 
The operators 2i , I2 and 73 then (up to an additive constant) take the form 

h = x(l-x-y)dl + (K 1 (y-l)-(l3 + K2 )x)d x , 
h = y(l-x-y)d% + (K 2 (x-l)-(j3 + Ki)y)d y , 
h = xy(d x - d y ) 2 + (n 2 x - K\y)(d x - d y ). 

Remark 3.2 The choice of variables entering the definition of G is dictated by the potential U , but 
the specific form needs to be determined. 

If we transform to case IX the second order parts of the two L operators are related in the usual 
way, but the first order parts take the form 

2ki + 1 + (2/3 - l)x 2 2k 2 + 1 + (2/3 - l)y 2 

OX H Oy, 

x y 

containing non-polynomial parts. To remove these we must set k% = k 2 = — i. The other coefficients 
are related by (3 9 = 2(3\ — 1. It should be noted that this gives the full parameter version of Case 
IX fromn a restricted version of Case I. 

3.2 Super-Integrability for Metrics II and III 

If we consider L = L b + U, h = H 2 + V 1: I 2 = E 2 + V 2 , then we find 

tt fc i , fc 2 , k 3 {y-l) T/ 16fci(l-y) 16fc 3 (l-y) 2 T/ x y 

U= 1 1 - , V\ = - , V 2 = -ik 2 4fc 3 -. (10) 

x y x z x ar y x 
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Remark 3.3 The commutator [7i, 7 2 ] is a third order operator, so cannot be written as a polynomial 
expression in 7, I\, 7 2 . However, the next commutators, [7i,/ 2 ]] o-nd [7 2 , [7i,/ 2 ]] can, giving us 
a quadratic algebra (see W), but this will not be considered here. 

Under the involution, these operators transform to 7 = 7& + U, I\ = H 2 + Vi, I3 = F 2 + V3, 
where 

- k x x k 2 k 3 x 2 - 16fcia; 16k 3 x 2 (y-1) 2 xy 

U = 7 77^ H ^7 7T7j V\ = yz r - 7- ttt, V3 = -4fc 2 4fc 3 



(y-1) 2 y (y-1) 3 ' (1-2/) (1-y) 2 ' ~ xy "(y~l) 2 ' 

If we consider L = L b + U, I 2 = E 2 + V 2 , I 3 = F 2 + V 3 , then we find 

u= k A v 2 = ^, F3 = -4fc 2 (y-D 2 ; 

y y xy 

and, furthermore, [H, L] = 0. These operators all commute with L and form the algebra: 

[77, I 2 ] = 47 2 , [77, I3] = -47 3 , [7 2 , J 3 ] = 16/77 - 77 3 + 2(8k 2 - 1)77. 
They also satisfy the constraint: 

7 2 7 3 + 73/2 = 327 2 - 16(4fc 2 + 1)7 - 4i7 2 7 + -77 4 + -(5 - 8k 2 )H 2 + 16k 2 (2k 2 - 1). 

2 

3.2.1 The Transformation to Case III 

All of these formula can be transformed to Case III. For example, this last algebra is satisfied by 

where 7;,, 7 1 , F and i7 arc those in Case III. The numerical factor follows from the fact that £2 n- 
— 16753, F 2 — j$Fa, i7 2 i-> 77 3 , which satisfy the same commutation relations, since such a 
scaling is just another isomorphism of the algebra. 



3.2.2 The Krall-Sheffer Gauge 

Starting with the operator 7 = Lb + U, with U defined by (|10j) . the gauge transformation 

7 1 y G~ l LG, with G = e a[ :»- 1 ^ x x h y c 

leads to fci = a(l — 2b), k 2 = \c(2c — 1), fc 3 = a 2 if we require the potential to reduce to a constant. 
The operator then reduces to the Krall-Shcffcr form with 

a=y, b=±(J3 + K 2 -l), c=~(2K2 + l). 

The operators 1\ and 7 2 then (up to an additive (and an overall multiplicative) constant) take the 
form 

h = x 2 d 2 x + ((/3 + k 2 )x + «i(l - y))d x , I 2 = xyd 2 + (my - n 2 x)d y . 
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Transforming to case III: the second order parts of the two L operators are related in the usual 
way, but the first order parts take the form 

containing non-polynomial parts. To remove these we must set k,\ = 0, ft 2 = — 5, which leaves us 
with the reduced operator with first order part given by (3xd x + (3yd y . 

To obtain the most general Krall-Sheffer operator we need to consider the operators j(HE+EH) 
and E 2 . We find the following result: 

T ki k 2 y k 3 (x + y 2 ) 
L — L b H 1 — H 

X X z X 6 

commutes with the two operators 

h = 1 -{HE + EH) *lg - M * ± ^ - hy(2 \ + ^ and J 2 =^ + ^ + ^!. 
4 x x z x A xx 

the gauge transformation 

L 1 y G~ X LG, with G = x ^ e H^+y 2 )/x 2 +c y / x 

leads to k\ = h(2b(4a — 3) — 2c 2 ), k 2 = —Abe, k 3 — —4b 2 if we require the potential to reduce to a 
constant. The operator then reduces to the Krall-Sheffer form with 

a=Y> b=^(J3 + K2-l), c=--(2K2 + l). 
The operators I\ and I 2 then take the form 

I x = 2x 2 d x d y + xyd 2 + {k 2 x - K\y)d x + ((3x + Kt)d y , I 2 = x 2 d 2 + (n 2 x - Kiy)d y . 

3.3 The Metrics IV, VI and VII 

These cases are separated in Cartesian coordinates, so the examples of super-integrable cases are 
very simple and well known. In the Krall-Sheffer gauge the operators L are just direct sums of 
1— dimensional ones, corresponding to products of classical polynomials, so we omit these examples 
here. 

3.4 Super-Integrability for Metrics V and VIII 

If we consider L = L b + U, h = H 2 + V x , I 2 = E 2 + V 2 , then we find 

U = k 1 y-\ 1 -, Vi = 5-, V 2 = kix . (11) 

<l ■ >> ■ >A < / • >> ■ 

The involution transforms to the case with L = L b + U , I\ = H 2 +V\, I3 — F 2 + V 3 . 
We may also consider the case L = Lb + U, I2 = E 2 + V2, I3 — F 2 + V 3 , which gives 

_ k 3 y _ k\y 2 2k2y 
U — kiy + k 2 vx + —j=, V 2 = hx + 2k 3 ^/x, V 3 = 1 

y2) X -y X 

which is invariant under the involution. 
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3.4.1 The Krall-Sheffer Gauge 

Starting with the case (jlip . the gauge transformation 

L^G- l LG, with G = e (ax+b ^' x x c 

leads to k\ = —a 2 , k 2 = b(l — 2c), ^3 = b 2 if we require the potential to reduce to a constant. The 
operator then reduces to the Krall-Sheffer form with 

(3 , Kl («2 - 1) 

d = — , = — , c= . 

2 2 2 

The operators I\ and -Z2 then (up to an additive (and an overall multiplicative) constant) take the 
form 

h = x 2 d 2 + (n 2 x - Kiy)d x , I 2 = xd 2 + (f3x + Ki)d v . 

Case VIII: Transforming between these cases restricts the parameters to the Laplace-Beltrami 
case. Here we need the choice L = L b + U, h = F 2 +V X , I 2 = \{E 2 - 2(HF + FH)) + V 2 , where 

U = ^y(y 2 - 2x) + ^(y 2 - x) + \(k\ - 20m), 

v x = ~(fi y + k 2 ) 2 , v 2 = ^(y 2 - x ) x + |( k 2 + m^y - m)). 

Using G = e i& K 2x+6K iy +6i3xy-3K 2 y 2 -2i3y 3 )/i2^ ^ggg g au g e transform (up to additive constants) to: 
L = ydl + 2d x d y + {(3x + K 1 )d x + {f3y + K2 )d y , 

h = d 2 + (fly + n 2 )d x , (12) 
h = (y 2 -x)d 2 x + 2yd x d y + d 2 + (K 1 y-K 2 x)d x + (l3x + K 1 )d y , 



4 The Polynomial Eigenfunctions 

By Definition 11.11 an admissible operator possesses a sequence of eigenvalues A„ (with A„ 7^ A m for 
m =/= n), such that, for each n > 0, there exist n + 1 linearly independent polynomial eigenfunctions 
of degree n, with eigenvalue A„. The degeneracy (for n > 0) of A„ stems from the super-integrability 
in a very explicit way. We show below how to use the operators I\ and I 2 to build sequences of 
polynomial eigenfunctions of L. 

We may choose a basis of monic polynomials P m ,n of the form ([2]) and organise them in a 
triangular array (see Figure [1]), with the degree N = m+ n constant on horizontal layers. Since the 
eigenvalue ([3]) depends only on the degree, it is constant on each horizontal level. 

Remark 4.1 Cases IV, VI and VII separate into sums of elementary 1— dimensional operators 
so their eigenfunctions are just products of classical orthogonal polynomials (respectively Laguerre- 
Laguerre, Laguerre-Hermite and Hermite-Hermite), so will not be discussed further. 

The Krall-Sheffer operators reduce to one dimensional operators on at least one of the edges of 
this triangle (usually both). The two exceptional cases are III, which only reduces on the left and 
VIII, which only reduces on the right. This means that there exist 3 point recurrence relations for 
the eigenfunctions P m ,o and/or Po. u . Below, we give explicit 3 level recurrence relations for each 
case and show how these reduce to the 3 point recurrence relations on the appropriate edges. The 
one dimensional reductions on the left and/or right edges also possess differential ladder operators. 
We also present the ladder operators in directions which are "parallel" to the edges for each case 
and show how these reduce to their one dimensional counterparts on appropriate edges. 
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What is the role of the commuting operators? Since 

Lf = Xf 1 
[/,£]= J 



L(If) = X(If), 



the function // is also an eigenfunction with the same eigenvalue. For our Krall-Sheffer operators, 
we have two such operators I± , I2 which take us respectively right and left on each horizontal level. 
Starting from Pn,o or Po,at, we use these operators to build the eigenfunctions P m , n , with N = m + n. 

r P °'° r 

/ ^ 
Pi,o ^ Pq.i 

, /_ \ 

Pifi <^ Pi,i <^ P0.2 

, / \ 

■P3,0 <^ Pm <^ -Pi. 2 ^ -Po.3 

/ \ 

P4.O <^ p3,l <^ p2,2 Pi, 3 <^ Po,4 

Figure 1: The triangular lattice of polynomials P m ,n with Po.o = 1- Horizontal arrows denote the 
action of 1\ (right) and I2 (left). The operators and r y represent the edge ladder operators. 



The 3 Level Recurrence and Ladder Operators. It is possible to explicitly determine the 
coefficients of the 3 level recurrence relations (|4]). This is a laborious calculation, so only the resulting 
recurrence relations are presented. The general form of Q constitutes a difference equation on 
2(m + 71+1) points, but our formulae contain at most 9 points, and each equation fits into the 
following 12 point relation: 

3 6 
Pm+l,n ^Pm,n + ^ ^ ^'-^ m-l-H,n+l-i ~t~ ^ ^ ^%Pyn— 3+i,n+2— i: 



i=0 i=0 

m+n 6 

Pm,n+1 yPm.n + ^ ^ CiPm+1— 1-H + ^ ^ ^iPm+2 — i,n— 3-H 7 

i=0 i=0 

Given the first of these, we operate with 

£ = ^4<9^ + 2Bd x d y + C&1 + Pd, + Ed y , 

taking into account that P m , n , is an eigenfunction, with eigenvalue A m +„, we find 

3 



(13) 



X m +n+lPm+l,n — X m -\- n (xP m n + S ' fli-Pm-l + i,n+l-i ) + (2.45, + 25^ + £>)P„ 

i=0 

6 

+ \n+n-l biPm-3+i,n+2-i- 
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Subtracting A m + n _i times the first of (|13[) , leads to a first order two level differential-difference 
relation: 

3 

(Am+n+l ^m|Ti-l)-fm+l,n {\n-\-n A m _|_ n _ l) (^m,Tl "f" ^ ^ OiPm — lH-i,ra+l — i) 

i=0 

+ (2Ad x + 2Bd y + D)P m . n . (14) 

The next step in the calculation is to use the horizontal shift operators I\ and I 2 to write all the 
polynomials of degree m + n (arising in (|14[1 ) in terms of the polynomial P TO , n . In principle, this 
could result in a fourth order differential operator, but only second and first order operators occur. 
This, together with a similar manipulation of the second of (|13[) . gives two differential operators 
which raise the degree of the polynomial eigenfunction: 

It is also possible to construct lowering operators, but these are usually considerably more compli- 
cated. 

These operators (usually) depend explicitly upon both m and n (not just through the combination 
m + n). The formulae for [L,P +i J and [L,R+ tV ] are non-standard but still imply that P_|_ jX P m;n 
and R+ t yP mtn are eigenfunctions of L, with eigenvalue A m + n +i. 

4.1 Krall-Sheffer Metric I 

Here we have 

L = (x 2 - x)dl + 2xyd x d y + (y 2 - y)d 2 + (fix + n{)d x + (J3y + K 2 )d y , 
h = x(l - x - y)d 2 x + (Ki(y - 1) - (/3 + K 2 )x)d X) 
h = y(l-x-y)d 2 + (K 2 (x-l)-(l3 + K 1 )y)d y , 
I 3 = xy(d x - d y ) 2 + (k 2 x - niy){d x - d y ). 
The polynomial eigenfunctions of L satisfy the 6 point difference equations: 

Pm+l,n %Pm,n 

+ (((/? + 2n - 2)(«i - 2m) - 2m(m + l))P m ,„ + 2n{n - k 2 - l)P m+ i,„_i) 

72W72JV-2 

H J" ( m C# + m + 2n- 2)(«i - m + l)(/3 + «i + m + 2n - l)P m -i n 

72AT-172JV-272JV-3 

+ (n(/3 + 2n - 3)(ki - 2m)(K 2 - n + 1) - 2m(?7i + l))P TOi „_ 1 

-n(n - 1)(k 2 - n.+ 1)(«2 -n+ 2)P m+ i jn _ 2 ), (16) 

and 

Pm,n+1 yP? 



1 



(((/3 + 2n - 2)(ki - 2m) - 2m(m + l))P m ,„ + 2m(m - /ti - l)P m _i, n+ i) 



72Af72iV-2 

2- ( n (/3 + 2m + n - 2)(re 2 -n+ l)(/3 + ft 2 + 2m + n- l)P m „_i 

72JV-l7 2W -272W-3 

-(m(/3 + 2m - 3)(k 2 - 2n)(«i - m+ 1) - 2ra(ra + l))P m _i in 

-m(m - 1)(ki - m + 1)(k x - m + 2)P m _ 2 , n+ i), (17) 
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where N = m + n, j n = ft + n. The 6 points of these difference equations are depicted by bullet 
points in Figure [2j Notice that when n = 1, the point (m + l,n — 2) is out of the Figure, but also 

o 

o o 
o o o 
o o o o 
o • • • o 
o o • • o o 
o o o • o o o 



Figure 2: The pattern for the difference equations for both P m +i,n and P m ,n+i for both cases I and 
II. 

that in formula (fT5|) the coefficient of P TO +i,n-2 vanishes when n = 1, so the formula is still valid. 
When n = 0, the three points (m + 1, n — 2), (m + 1, n — 1) and (m, n — 1) fall outside the Figure, 
but again, the corresponding coefficients vanish, leaving us with the 3 point recurrence relation on 
the left edge: 

p ( (/3-2)(« 1 -2m)-2m(m+l) \ 

ro+1 »° ~ V (/3 + 2m)(/3 + 2m-2) 



m(/3 + m — 2)(ki — m + l)(/3 + «i + m — 1 
(/3 + 2m - l)(/3 + 2m - 2) 2 (/? + 2m - 3) 



Pfn-lfi- (18) 



Similarly, as we approach the right edge, the formula (flT|) consistently reduces to the 3 point recur- 
rence relation on the right edge: 

_ ( (P-2)(K 2 -2n)-2n(n + l) \ 
P °'" +1 " [ y+ (/3 + 2n)(/3 + 2n-2) J P °-" 

rtQg + n - 2){k 2 - n + l)(/3 + k 2 + n - 1) \ 
(/3 + 2n-l)(/3 + 2n-2) 2 (/3 + 2m-3) / 1 J 

The left and right edge polynomials satisfy respective Jacobi equations: 

L [x) P mfi = (x{x - l)d 2 x + (fix + Ki)d x )P mfi = m(m- l + /3)F m , , 
LMp , n = (y(y-l)dl + (l3y + K 2 )d y )P 0tn = n(n - 1 + /3)P ,„. 

To proceed with these second order difference equations, we expect to need the three lowest polyno- 
mials: 

fo,o = l, Pi,o = x+j, P 0A = y+j, (20) 

which are calculated directly as monic, degree one polynomial cigcnfunctions. However, when m = 0, 
(|18|1 is consistent and can be used to construct Pi t o (and similarly for Po,i)- We then use the difference 
equations (fl~6f and ([T7|) to calculate the degree 2, degree 3, etc polynomials successively. 

Alternatively, we may calculate the polynomials P m ,o and Po,n; using (JTSJ) and (|19p and then 
use the commuting operators to calculate the remaining polynomials on each level. We know that, 
since the operators I\ and I 2 commute with L, the functions IkPm,n,k = 1,2, are also polynomial 
cigcnfunctions of the same degree. By looking at the leading order (monic) part, we obtain: 

hP m ,n + m(fi + k 2 + m - l)P m ,n = m(«i -m+ l)P m _i jn+ i, 
hP m ,n + n(/3 + ki + n - l)P m ,n = n{n 2 - n + l)P m+ i n _i, 
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We see that, starting with P m ,o, I\ moves us to the right, until we reach Po,n, which satisfies 
IiPa,n = 0. I 2 similarly moves us from right to left. 

Finally, we may construct the operators (|15[) . which are given by 

R+,x = ^ - 2N){ l +2N _ X) ((P + N- 1)(08 + 2N)x + Kl -N) + (/3 + 2N)x(x - l)d x 

+ ((f3 + 2N)xy + {fi+ m)y + k 2 (1 - x))d v + y(x + y- 1)^), 
1 

(P + 2N)(0 + 2N- 1)' 



R+.y = la , , nAT ,M P + N-l)((P + 2N)y + K2 -N) + (P + 2N)y(y-l)d y 



+ (((3 + 2N)xy +(j3 + k 2 )x + - y))d x +x(x + y- 1)%). 

The formulae interchange under the transformation x <-> y,m <-> n 2 . These operators satisfy the 
commutation relations 

(2x-l)(L-X N ) 

[-^j-K+.xJ - /g + 27V - 1 r l A W+l ~ 'V/Vj-H.+.x) 

(2y-l)(L-A JV ) 
[-k> - — /3 + 27V - 1 ( A ^+i - *N) H +,y 

Starting from Po,n, is used to build the functions P m ^ n ,m > 1, with n fixed. When n = 0, 
reduces to the 1 dimensional ladder operator 

r+. x = ^ + 2 m)(l+2m-l) ^ + ™ ~ + 2m)x + Kl ~ ^ + (/? + 2m ^ X ~ 

Similar statements can be made regarding P+,y. 

This particular case is symmetric in (x,y) and (/ci,K2), which is reflected in all the formulae. 

4.2 Krall-Sheffer Metric II 

Here we have 

L = x 2 dl + 2xyd x d y + (y 2 -y)dl + (f3x + Kl )d x + (f3y + K2 )d y , 
h = x 2 dl + ((/3 + K 2 )x + Kl (l-y))d x , 
I 2 = xyd 2 + (my - n 2 x)d y . 

The polynomial eigcnfunctions of L satisfy the 6 point difference equations: 

Pm+l,n = %Pm,n 

+ (m(P + 2n- 2)P m . n + 2n(n - k 2 - l)P m+ i,„_i) 

72W72AT-2 

H 2 (mn\(l3 + m + 2n- 2)P TO _x n 

72W-172W-272W-3 

+(nm(/3 + 2n~- 3)(« 2 - n + l))P m ,„_i 

-n(n - 1)(k 2 - n + l)(«a - n + 2)P m+ i j „_ 2 ), (21) 
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and 

(((/3 + 2m - 2)(«2 - 2n) - 2n(n + l))P ro ,„ - 2mKiP m -i,n+i) 



J2NJ2N-2 

H 5- (n(/3 + 2m + n - 2)(k 2 - n + l)(/3 + K 2 + 2m + n - l))P m „_i 

72AT-17 2 at_ 2 72AT-3 

+(m/ci(/3 + 2m - 3)(k 2 - 2n) - 2n(n + l))P m _i,„ - m(m - l)^P m _ 2 „ +1 ), (22) 

where N = m + n, 7 n = /3 + n. The 6 points of these difference equations are again depicted by the 
bullet points in Figure [5] Again, as we approach the left edge, the formula (|2"Tj) consistently reduces 
to the 3 point recurrence relation: 

Pm+1 >° ~ [ X+ (/3 + 2m)(/3 + 2m-2)J Pm '° 

mnf((3 + m — 2) 

Pm-i,o- (23) 



03 + 2m - 1)(J3 + 2m- 2) 2 (B + 2m - 3) 

Similarly, as we approach the right edge, the formula (|22[) consistently reduces to the 3 point recur- 
rence relation: 



, (B-2)(n 2 -2n)-2n(n + l) \ 

y+ (/3 + 2n)(/3 + 2n-2) ) P °<™ 

n(8 + n- 2)(k 2 -n+ 1)00 + k 2 + n - 1) 



08 + 2n - 1)09 + 2n - 2) 2 (/3 + 2n - 3) 
The left and right edge polynomials satisfy respectively: 

L^P mfi = {x 2 dl + (J3x + Kl )d x )P m v = m{m-l+P)P mfi , 



P m ,n-l. (24) 



(25) 



L^P 0in = (y(y-l)d 2 + (By + K 2 )d v )P , n = n(n - 1 + /3)P ,„. 

Again, we have the polynomials (f20|) . which are the seeds for the difference equations (|21|) and (|22|) 
when calculating the degree 2, degree 3, etc polynomials. 

Alternatively, we may calculate the polynomials P m ,o and Po.n, using (f23| and (|24|) and then use 
the commuting operators to calculate the remaining polynomials on each level. By looking at the 
leading order (monic) part of ifcP m , n , we obtain: 

hPm,n - m(J3 + k 2 + m - l)P m ,„ = -m/«iP OT _i in+ i, 

hP m ,n - nniP m . n = -n(« 2 - 7i + l)P m+ i, n _i, 

We see that, starting with P m ,Oi h moves us to the right, until we reach Po. n , which satisfies 
Ji-Po,n = 0. 1% similarly moves us from right to left. 

Finally, we may construct the operators (|15[) , which are given by 

R+,X = 7g + 2N){ p + 2N-l) i{l3 + N ~ lK(/? + 2N)X + ^ + ^ + 2N > 2 ^ 

+ ((/3 + 2N)xy + my - n 2 x)d y + xyd 2 ), 

R+ v = i w T ((B + N-l)((8 + 2N)y + K 2 -N) 

+ ' u [8 + 2N)(8 + 2N - )W >y 2 ' 

+ ((8 + 2N)xy + (B + k 2 )x + «i(l - y))d x + (8 + 2N)y(y - l)8 y + x 2 d 2 x ). 



16 



The formulae no longer interchange under the transformation x f-> y, K\ <-> k 2 - These operators 
satisfy the commutation relations 

[L,R+, X ] = p\ r 2N -1 + ^ JV+1 ~ ^ JV )-^+. a ' 

(2y-l)(L-Ajy) 
[-t>,it +jJ/ J - /g + 2iV - 1 (. A W+i - AN)-K+,y 

Starting from Po,n- P+.z is used to build the functions P m , n ,m > 1, with n fixed. When n = 0, 
R+ lX reduces to the 1 dimensional ladder operator 

r+ - x = (f3 + 2m)(J$ + 2m-l) {{P + ™ ~ + 2m)x + Kl) + {P + 2m ^ d ^- ^ 
On the right edge, P+, y reduces to 

T+ - y = (/3 + 2m)J+2m-l) ((/3 + " ~ + 2n)y + ^ n) + (fi + 2n)viV ~ W 

4.3 Krall-Sheffer Metric III 

Here we have 

L = x 2 dl + 2xyd x d y + (y 2 +x)dl + (f3x + ti 1 )d x + {Py + K 2 )dy, 
h = 2x 2 d x d y + xyd 2 + {k 2 x - K\.y)d x + {J3x + ni)d y , 
h = x 2 d 2 + (k 2 x - K\y)d y . 

The polynomial eigenfunctions of L satisfy the 9 point difference equations: 

P m +i.n = xP m ^ n H (ki(P + 2n - 2)P m . n - 2nn 2 P m +i,n-i - 2n(n - l)P m+2)n _ 2 ) 

72AT72AT-2 
1 



— {mKiiP + m + 2n — 2)P m _i i „ + nnin 2 ((3 + 2n- 3)P m ,„_i 



72JV-l72JV-272iV-3 

+n(n - 1)(kxC9 + 2n - 4) - K 2 )P m+ i.„- 2 - 2n(n - l)(n - 2)K 2 P m+2 .„_ 3 

-n(n-l)(n-2)(n-3)P m+3 ,„_4), (27) 



and 

Pn.n+l yPm,n 

+ M/3 + 2m - 2)P TOi „ + 2(n(/3 + 2m + n- l)P m+ i,„-i - m/siP m -i,„+i) 

72iV72iV-2 

+ J 1 (n(n - l)(n - 2)(2/3 + 4m + 3n - 3)P m+2 ,„- 3 

72jV-l7 2 AT_ 2 72N-3 

+n(n - l)(2/3 + 6m + 4n - 5)/t 2 Pm+i,n-2 + n({0 + 2m + n- 2){k\ - m(0 + 3n - 3)) 
+Ki(l - n 2 ))P m , n -i + mK!K 2 (/3 + 2m - 3)P m _i >n - m(m - l)K?P m -2,n+i, (28) 

where N = m + ?i, 7« = /? + n. The 9 points of these difference equations are depicted by the bullet 
points in Figure |3l 
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Figure 3: The pattern for the difference equations for both P„ l+ i „ and P mj „+i for case III. 



Note that in this case, the operator L does not reduce to a 1 dimensional operator on the right 
edge and, similarly, the difference equations cannot be reduced here. However on the left edge, the 
operator L reduces to L^ x > of ([25]) and the difference equation ([27]) reduces to (j23)) . Nevertheless, 
on the right edge, the difference equation ((28]) still makes sense, so can be used to construct the (2 
variable) polynomials on the right edge. 

Starting from the polynomials P m ,o we may use the commuting operators to calculate the re- 
maining polynomials on each level. By looking at the leading order (monic) part of IkPm,m we 
obtain: 

hP m ,n - mn 2 P m ,n = -m,KiP m -i,n+i+n(P + 2m + n-l)P m +i, n -i, 

hPm,n + riKlPm,n = nK 2 P m+ i y „-l + n(n - l)P m +2,n-2, 

which are each 3-point differential, recurrence relations on each level. At the beginning, when 
n = 0,m = N, the first gives a 2-point relation, giving Pn-i ,i- When m = 0, n = N, we have 
IiPq,n = N(f3 + N — 1)Pi ; jv-i. Having constructed the polynomials P m ,n, m + n = N in this way, 
we may operate with I2 to move to the left. In this case /2-P/V.0 = 0. 
Finally, we may construct the operators (|15[) . which are given by 

R+,x = (7^-1(72^2; + «i) + j2NX 2 d x + (-f 2 Nxy + K\y - n 2 x)d y - x 2 dl), 

72NJ2N-1 

R+.y = (jN-l (72NV + «a) + {liNxy + fax - Kiy)d x 

l2N"f2N-l 

+ (l2Ny 2 + Zjnx + Ki)d y + 2x 2 d x d y + xyd 2 ). 

The formulae no longer interchange under the transformation x y,Ki <-> K2- These operators 
satisfy the commutation relations 

[L,R+ X ] = — — + (Ajv+i — \n)R+.x, 

72JV-1 

[L,R + , y ] = — + (Aw+i — Xn)R+.v 

12N-1 

On the left edge, R+ :X reduces to the 1 dimensional operator (j26)) . which generates the polynomials 
Pm,o- On the right edge, R+. y does not reduce to a 1 dimensional operator, but still generates the 
polynomials Po,n- 
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4.4 Krall-Sheffer Metric V 

Here we have 



L = 2xd x d y +ydy + (0x + Ki)d x + (/3y + K 2 )d y , 



v 

„2a2 



h = x d x + (k 2 x - Ki_y)d x , 
h = xdy + (/3x + Ki)d y . 

The polynomial eigenfunctions of L satisfy the 5 and 4 point difference equations: 



Pm+l,n — [ X + — J P m . n H 7rP m +l,n-l — ~7P>(( n ~ l)-fm+X,n-2 + ^l-Pfre.n-l)) (29) 

P / P P 



and 



K2 + 2(m + n)\ 1 , 

■Pm,n+i = ( J/H o ) %,n - -^{n(K2 + 2m + n - l)P TO) „_i + mKiP m _i, n J. (30) 



The points of these difference equations are depicted by the bullet points in Figure [4] 

o o 

o o o o 

o o o o o o 

oooo oooo 

o • • o o o • • o o 

o o • • o o o o • o o o 

o o o • o o o o o o • o o o 



Figure 4: The pattern for the difference equations for (a) P m +i,n\ (b) P m ,n+i- 
On the left edge (f2"9")l reduces to 

( K\\ ( K\ 

Pm+1,0 = I X + — I P m ,0 => Pmfi = I X + — 

On the right edge, this reduces to (f30| reduces to 



Po,n+l = (y-\ a ) Po,n - -^(n{n 2 + n - l)P , n _i), 



which is the three point recurrence relation for the Laguerre polynomials. 
Notice that the operator L reduces to 

(fix + = mf3P m .Q, 

on the left edge, giving the simple solution above for P m ,o- On the right edge, the operator reduces 
to that of Laguerre. 

Starting from the polynomials P m ,o we may use the commuting operators to calculate the re- 
maining polynomials on each level. By looking at the leading order (monic) part of IkPm,n, we 
obtain: 

hPm,n — rn(K2 + m— l)P m ,n = ~ WKlP m _l )n +l, I 2 Pm,n + HKlPmn = nf3P m +l,n-l-i 
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with I\Po, n = and hPmfi = 0. 

Finally, we may construct the operators (|15[) , which are given by 

R+, x = {xd 2 y + (2px + Kl )d y + P{px + , 

R+, y = p (xd x + yd y + Py + m + n + k 2 ). 

These operators satisfy the commutation relations 



[L, R+, x ] — /3R+,x> [L, R+,y] — 1~ 0R+,v 



(L _X N ) 

P 

On the left edge, R+ lX reduces to the multiplicative operator 

r +,x — X -t — , 

giving the above simple formula for P m .o- On the right edge, R+, y reduces to 

r+, y = -b (ydy + py + n + k 2 ). 

4.5 Krall-Sheffer Metric VIII 

Here we have 

L = ydl + 2d x d y + {px + K 1 )d x + {py + K 2 )d y , 
h = d 2 x + (/3y + K 2 )d x> 

h = {y 2 - x)dl + 2yd x d y + d% + (my - n 2 x)d x + {fix + Ki)d v . 
The polynomial cigenfunctions of L satisfy the 5 and 3 point difference equations: 



Pm+l.n — [ X + — ) P m n + — (2mP m -l n +l + TlP m ,n-l) T^ - Pm-1. 

P J P P 



and 



n 2 \ TYl 
P?n,n-{-l — I y ~7T 1 Pra,n "T ~T P m — 1 , n • 



P ) """ P~ 

The points of these difference equations are depicted by the bullet points in Figure El 

o o 

o o o o 

o o o o o o 

oooo oooo 

o o • • o o o • o o 

o o o • • o o o • o o o 

o o o • o o o o o o • o o o 



Figure 5: The pattern for the difference equations for (a) P m +i,n, (b) P m ,n-\ 
On the right edge, (|32)) reduces to 



K2\ I K 2 

Po,n+l = [y+-g-J "0,n => Po,n = I V + -p 
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On the left edge (f3"Tj) reduces to 

/ Ki\ 2m mK2 

Pm+1,0 = I X + — P m o H ^- im— 1,1 T^T" °ro-l,0) 

\ P J P P 

which is no longer a 3 point recursion on the edge, since the left edge polynomials are not functions 
of just x. 

Starting from the polynomials Pq_ u we may use the commuting operators to calculate the re- 
maining polynomials on each level. By looking at the leading order (monic) part of IkPm,m we 
obtain: 

hPm,n = TnpPm-\,n+\, 

+ m(m — l)P m _2,n+2) 

with I\Po,n = 0. 

Finally, we may construct the operators (fT5|). which are given by 



K 2 , 1 



R+, v = y + -j + p d *- 

These operators satisfy the commutation relations 

[L, R+,x] = f3R+,x, [L, R+, y ] = (3R+,y 
On the right edge, R+, y reduces to the multiplicative operator 

r+,y =V+ -j, 

giving the above simple formula for Po, n - On the left edge, R+. x does not reduce to an operator in 
just x, but still generates the functions P m ,o from Pq^q = 1. 

4.6 Krall-Sheffer Metric IX 

In this case it is convenient to consider the 4 operators, which commute with L: 
L = (x 2 - l)fig + 2xyd x d y + (y 2 - \)d 2 y + (3xd x + f3yd y 

h = {l-x 2 -y 2 )dl + {l-p)xd x , I 2 = (l-x 2 -y 2 )d 2 + (l-(3)yd y , (33) 
h = xd y - yd x , h = 2(1 - x 2 - y 2 )d x d y + (1 - /3)(xd y + yd*). 

Under the commutator, these form a closed, quadratic algebra (see [10]). which we don't need here. 
They are not, of course, independent. They satisfy the quadratic relations: 

h + h + if + L = 0, 2{hh + hh) - W 2 - 4/3 - l)(h + I 2 ) =(fi- - 5)L + I 2 . 

These also satisfy 

hP m ,n = (2 - P - m)mP mm - m(m - l)P m -2,n+2, 

hPra.n = (2 - /J - Tl)nP m , n - fl(n - l)P m+2 ,n-2, 
4fm,n = n Pm+l,n-l ~ m Pm-l,n+li 

hP m . n = (1 - P ~ 2m)nP m+ i ! „_i +(!-£- 2n)mP m _i <n+1 , 
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so we use -Z3 and I4 as our horizontal ladders. 

The polynomial eigenfunctions of L satisfy the 4 point difference equations: 

P m +i,n = xP m ,n H (n(n - l)P m+ i,„_2 - m(/3 + m + 2n - 2)P m _i n ), 

72JV-172JV-3 

and 

Pm,n+i = yP m ,n H (m(m - l)P TO _2, n +i - n(J3 + 2m + n — 2)P m . n _i), (34) 

72IV-172JV-3 

where N = m + H, 7« = /? + n.. The points of these difference equations arc depicted by the bullet 
points in Figure |6l On the left edge, the first reduces to 

Pot+i,o = xPm,o m(j8 + m - 2)P m _i. , 

72m-l72m-3 

which is the 3 point recurrence relation for Gegenbauer's equation 

/1 2\^ 2 ^ 5 "i,o a dP m o , , 

[1 - x )— — px— - — = m(m + p- l)P ro ,o- 

dx z ax 

The second relation (|34|) similarly reduces on the right edge. 

o o 

00 00 

000 000 

0000 0000 

o • o • o o • o • o 

ooo»oo oo»ooo 
ooo»ooo ooo»ooo 

Fi gure 6: The pattern for the difference equations for (a) Pm+i,m (b) Pm,n+i- 
Finally, we may construct the operators (fT5|). which are given by 

R+,x = —?—(xyd y + (x 2 -l)d x + (P + N-l)x), 

72JV-1 

R+, y = -^—(xyd x + (y 2 -l)dy + (/3 + N-l)y). 

72JV-1 

These operators satisfy the commutation relations 

[L,R+ >X ] = — — + (Ajv+i — ^n)R+,xi 

72JV-1 

[L,R+. y ] = — — — + (Ajv+i — Ajv)P+,j/- 

72W-1 

5 Generating Functions 

Generating functions are known for the special cases IV, VI and VII, since these are just products of 
the known classical generating functions. The generating function for Case IX is just a direct gener- 
alisation of that for Gegenbauer polynomials (in 1 dimension). A multi-dimensional generalisation 
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can be found in |3J. We believe that the generating functions given below for Cases V and VIII are 
new. 

Generating functions are important for a number of reasons and can be considered as fundamen- 
tal. Suppose a given family of polynomials P m , n are coefficients in a generating function expansion 

G{x^y^S^t^j y Amn-^mn \ T; 

* — ' mini 

m,n— 

where the additional coefficients A m _ n arc included to allow our polynomials to be monic. Then, by 
considering ^ and we can construct explicit forms of difference equations, relating polynomials 
of degrees m + n + l,m + n and m + n— 1. By considering ^ and 4^, we can construct explicit forms 
of differential ladder operators. Manipulating these, we can construct the eigenvalue problem for the 
polynomials, thus proving that they really are eigenfunctions of a particular differential operator. 



5.1 Cases I and IX: Generalised Gegenbauer Polynomials 



Since the operator L of (|33|). for Case IX, is symmetric in x and y it is easy to generalise the standard 
Gegenbauer generating function to 

G(x,y,s,t) = (l-2 S x-2ty + S 2 + t 2 Y 1 -^/ 2 = £ ^ 



m.n—0 



Notice that this generating function has a pair of symmetries (s, x) 1— > (— s, —x) and (i, y) 1— > (— t, —y), 
which must be reflected in each term P mn s m t n . When m is even (odd), P m ,n is an even (odd) 
function of x (and similarly for n and y). This means there will be an even polynomial of x, 
P2m.m such that ?2m+i,n = xp2m.n (with a similar statement for P m .2n+i)- Indeed, the first few 
polynomials are 

2 1 2 1 

P fl = 1, Pi,o=x, P 0ll =y, P 2 ,o = x - — — - , P M = xy, P ,2 = 2/ 



1 + /3' x,± u "~ y 1 + /?' 



P 3 , = a; ( x 2 - — Z — ] , P 2 ,i = y f .t 2 



A,2 = a; ( y 2 - ^-ru J . p o, 3 = y(y 2 - 



3 + pJ' uw " V 3 + P 

In Figure [7] below, polynomials which are even in both x and y are depicted by black dots. They 
occur only when N = m+n is even. This configuration is unimportant in the context of case IX, but 

• 

o o 
• o • 

0000 
• o • o • 
000000 
• o • o • o • 



Figure 7: The pattern of polynomials for case IX. The black dots are when both m and n are even, 
when the transformation (xg, j/g) = {y/x\, ^Jyi) (see Section |2~T|) is used, the black dots represent the 
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polynomial eigenfunctions in Case I (with the restriction K\ = K2 = — 5). Generally, P m , n (of Case 
IX) maps to Psl^h (of Case I), giving us "half- integer" indices, corresponding to non-polynomial 

eigenfunctions. For example, P2,\ Pi, 1 = \fy (x — 2 (p+i 

In Section [4761 we chose to use 13,14 to take us right and left on horizontal levels, since they 
just related "nearest neighbours". The operators jump two steps at a time, connecting (in 

particular) the black dots on the horizontal levels. These transform to the Ji, I2 given in Section |4~T1 
(with Ki = K2 = —\ and the replacement /3g >->■ 2/?i — 1). With these values of «;,, we have 

hP m ,n + m + m - |^ P m . n = m Q - P m _i, n+ i, 

so, when m = i, the right hand side vanishes and Pi n is just an eigenfunction of I\. This means 
that starting on one of the half-integer functions, I± moves us to the right until we reach the last 
one on the given level. The operator I2 acts in a similar way. The operators 13,14 can also be 
transformed to the reduced Case I, but these operators no longer have polynomial coefficients. 

5.2 Cases V and VIII 

For Case V, the generating function is 

f p 2 sx + ( Kl s + f3ty)((3 -t) \ ( ty K2 ~ s m t n 

G(x,y,S,t) = exp( — ) (1--) = P ™,n{x,y) — j-f. 



m.n— 



When £ = this reduces to Gi(x,s) = exp + and generates the polynomials P mj o 

\ m 

x+ ^) j located on the left edge. With s = 0, the generating function reduces to 



Gr {y,t) = e^\ 7 pL-\ (l- 1 



[fi-t)j V P 

and generates Laguerre polynomials P n on the right edge. In the interior of the triangle the 
polynomials depend on both x and y in a non-trivial way. For example 

p " = ( i+ j) (" + jH(* + § 

By differentiating G(x, y, s, t) with respect to s we find 

r)C 

((3-t) 2 ^ = (p 2 x + Kl ((3-t))G, 

OS 

which leads to the 3 level recurrence relation 

f K\\ 2n n 

"m+l,n = I X + — I P m ,n + —pPm+l,n-l — ~p2^ n ~ l)Pn+l,n-2 + KlPro,n-l)i 



which we already saw in Section 14741 When n = 0, this reduces to P m +i,o = ( x + ) P m .o, which 
gives the above P m o- 
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Differentiating G(x,y, s,t) with respect to t leads to a more complicated formula, but can be 
combined with the s— derivative to give 



Bf F)C 

9 - *)* «T - 2s (P - t )jT = W*V + K 2(/3 -t)- Kl s)G, 
at as 



which leads to the other 3 level recurrence relation 



K2 + 2 (to + n) \ 1 

fm,n+l = I 2/ H o I -Pm,« ~ 03 («(«2 + 2to + 71- l)P m ,„-l + TOKi_P r „_i in ), 



also given in Section 14.41 

For Case VIII, the generating function is 



G(x,y,s,t) = exp 



4s 3 + 3(2/3y + n 2 )s 2 + 6fist + Qp{tfx + «]> + 0% + K 2 )t) 



^2 p m A x ,y) 



6/3 2 

s m t n 



m\n\ 

m.n— 



Similar manipulations as described above lead to the difference equations (|3T|) and (|32j) . 

The transformation between the Laplacc-Bcltrami operators of Cases V and VIII (see Section 
12. 4[) do not extend to the Krall Sheffer operators in this case, so the polynomials generated for either 
case cannot be transformed to eigenfunctions of the other in this case. 



6 Conclusions 

In this paper we presented a number of recursive algebraic structures which allow the explicit 
construction of polynomial eigenfunctions of the Krall-Sheffer differential operators. By arranging 
these polynomials in the triangular array of Figure HJ we showed that the 2— dimensional Krall- 
Sheffer eigenvalue problem reduced to a standard 1— dimensional problem on at least one of the edges, 
where we could use the standard 3— point recurrence relation to build polynomial eigenfunctions of 
one variable. We then used the commuting operators (coming from super-integrability) to build the 
2— dimensional polynomials on the horizontal level of the triangle. 

We then gave explicit formulae for the 3— level recurrence relations for each of the (non-trivial) 
Krall-Sheffer cases. These reduce to the corresponding 3— point relations on appropriate edges. 
The 3— level recurrence relations give a very effective way of building the polynomials. We also 
gave explicit formulae for the differential raising operators, which shift us to the next eigenvalue 
in directions parallel to the edges of the triangular array. Lowering operators also exist (see [10]) 
and sometimes have simple formulae. However, for some of the cases, these operators are very 
complicated, so we didn't include any of them here. 

Finally, we discussed the existence of generating functions for the Krall-Sheffer polynomials. 
We presented new generating functions for cases V and VIII. It was shown how the well known 
generalised Gegenbauer generating function of Case IX can be used to build eigenfunctions of Case 
I with special parameter values. It would be interesting to generalise this to a generating function 
for the full Case I. Generating functions for Cases II and III are also currently unknown. 

The discussion of Sections [2] and [3] give us a direction to follow if we wish to build Krall-Sheffer 
like operators in higher dimensions. Some calculations in this direction can be found in |10j . but it 
is an important subject for future research. 
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